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Abstract: We show that the one dimensional unitary matrix model with poten- 
tial of the form all + bU^ + h.c. is integrable. By reduction to the dynamics of 
the eigenvalues, we establish the integrability of a system of particles in one space 
dimension in an external potential of the form acos(a: + a) + 6cos(2x + j3) and in- 
teracting through two-body potentials of the inverse sine square type. This system 
constitutes a generalization of the Sutherland model in the presence of external 
potentials. The positive-definite matrix model, obtained by analytic continuation, 
is also integrable, which leads to the integrability of a system of particles in hyper- 
bolic potentials interacting through two-body potentials of the inverse hypebolic 
sine square type. 
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In one space dimension an integrable class of systems is known, involving particles 
coupled through two-body potentials of a particular form. The generic type of these 
potentials is of the inverse square form. Calogero [1] first solved the three-body problem in 
the quantum case for inverse square interactions and quadratic external potential. Later 
on, the full iV-body problem was solved and shown to be integrable in both the classical 
and the quantum case [2] and further to be related to Lie algebras [3] . Sutherland [4] solved 
the problem with inverse sine square interactions (and no external potentials), which can 
be thought as the inverse square potential rendered periodic on the circle. Eventually, 
it was realized that the system of particles with two-body potentials of the Weierstrass 
function type is integrable [5,6]. Again, these potentials can be thought as the inverse 
square potential rendered periodic on a complex torus. For a review of these systems and 
a comprehensive list of references see [6] . 

An interesting feature of the above systems is that, at least some of them, admit a 
matrix formulation [6,7]. Specifically, the inverse square potential arises out of a hermitian 
matrix model, and the inverse sine square potential arises from a unitary matrix model. 
The matrix formulation is in many respects a better framework to study these systems. 
Such hermitian matrix models have been studied in physics in the context of large-iV 
expansions [8-10] and, recently, non-perturbative two-dimensional gravity [11]. Although 
the one dimensional (c = 1) unitary model has not been studied in this context, discrete 
(c < 1) models have been considered [12]. Further, the inverse square potential was shown 
to be of relevance to fractional statistics [13] and anyon physics [14]. It becomes, therefore, 
of interest to look for integrable systems with more general potentials. 

In a previous paper [15] we achieved such a generalization for the Calogero system; 
specifically, it was shown that the system of particles with inverse square potential in- 
teractions remains integrable at the presence of external potentials which are a gen- 
eral quartic polynomial in the coordinate. In this paper, we obtain a generalization of 
the Sutherland model; that is, we show that the system of particles with inverse sine 
square interactions remains integrable at the presence of external potentials of the form 
acos(x + a) + 6cos(2x + (3). An appropriate scaling limit of this system, then, is shown 
to reproduce the previous quartic system. In addition, the system with all trigonometric 
functions in the potentials replaced by their hyperbolic counterparts (i.e., the inverse sinh 
square system with hyperbolic external potentials) is also integrable. 

2 



The fact that integrability seems to work only for the above type of external potentials 
is somewhat puzzling. It is known, for instance, through a collective field description of 
the hermitian matrix model, that the large- N inverse square system is integrable for any 
external potential [16]. (Although there is no corresponding result for the unitary model 
we have little doubt that the same is true there.) This could be, though, just another 
special property of the large- N system. At any rate, the method of this paper seems to 
work only for the above-mentioned potentials. 

The system to be considered is a unitary matrix model in one time dimension with 
lagrangian 

C = -|tr([/" 1 [>) 2 - trW(U) (1) 

where U is a unitary N x N matrix depending on time t, and over dot denotes time 
derivative. The potential W(U) must be hermitian for consistency. The equations of 
motion from (1) read 

A {U~ l U) - W'U = (2) 



dt 

Due to the invariance of (1) under time independent unitary transformations of U, there 
is a conserved traceless matrix, namely 

[U.U-^^iP (3) 

as can explicitly be checked using (2). P is the generator of unitary transformations of 
U and constitutes a kind of conserved "angular momentum" in the (curved) space U(N). 
The eigenvalues of U, on the other hand, written as e lXn , n= 1, . . . N, can be thought as 
coordinates of N particles on the circle. Take, now, the particular case where all but one 
of the eigenvalues of P are equal, that is 

P = a(N\u >< u\ - 1) (4) 

where \u > is a constant A^-dimensional unit vector. This P is naturally obtained by 
gauging the U(N) invariance and coupling the system to fermions [17]. Then in can be 
shown with a method analogous to [17] that the eigenvalues of U satisfy the equations of 
motion 

l , S. — "V t* LOS rj 

x n = -V x n ) + ]T — 37FS 5 
^ 4 sin J Xn <> Xm 



where the potential of the particles is defined 



V(x) = W{e lx ) (6) 

and prime denotes x-derivative. These are exactly the equations of motion of particles 
of unit mass moving in the external potential V and interacting through the two-body 
potential 

V *W = TT^2 ( 7 ) 

In the abscence of external potentials, the above system is the Sutherland system which is 
known to be integrable. It is the purpose of this paper to show that the system remains 
integrable in the presence of external potentials V(x) of the form 

V(x) = a\ cos x + C12 sin x + b\ cos 2x + 62 sin % x (8) 

The above potential can be thought as a generalization on the circle of the quartic potential 
on the line. Indeed, the quartic potential is in fact a special case of the above potential and 
can be obtained from it in the limit of infinite radius of the circle. To see this, introduce 
explicitly the radius R by rescaling 

In terms of the new variables the potential becomes 

v(x) - jpvQ (io) 

Then, by choosing the coefficients of the potential to scale as 

01 = --R A a - 8R 6 c , a 2 = 2R 5 b 

hi = 2R 6 c + ^R 4 a , b 2 = -R 5 b (11) 
we see that upon taking the limit R — > 00 we recover the potential 

V(x) = ax 2 + bx 3 + cx 4 (12) 

It should be clear that this is essentially the only scaling that leads to a finite potential 
at the large R limit. Therefore, the integrability of the present model also contains the 
integrability of the quartic system as a limiting special case. 
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We shall then consider the unitary matrix problem with potential 

W(U) = ^U+ ^-V 1 + |[/ 2 + |V 2 , A = ai - ia 2 , B = b x - ib 2 (13) 

Through a rotation of the circle, that is, through a redefinition of U of the form U — > e^U, 
we can always redefine the phase of B. So, for simplicity we will choose B real. We shall 
define the "left-hamiltonian matrix" H as 

H = -\{U- 1 U) 2 + W{U) (14) 

whose trace gives the energy of the system. Similarly, we can define the "right-hamiltonian 
matrix" 

H = -\{UU~ 1 ) 2 + W(U) = UHU^ 1 (15) 

Although the trace of the above matrices is a conserved quantity, neither of them is con- 
served as a matrix, nor are their eigenvalues. Consider, now, the matrices 

L = H + -[U- 1 U,U + U- 1 } , M = -\U- l U + -(U - U- 1 ) (16) 

Then, upon using the equations of motion it is possible to show that 

L+[L,M] = (17) 

provided that we choose k as 

k 2 = B (18) 

Therefore, the matrices L and M constitute a Lax pair [18] and the eigenvalues of L are 
conserved. Equivalently, the traces 

In = trL n (19) 

are constants of the motion. Restoring, now, an arbitrary phase to B, through a phase 
redefinition of U, we get that the Lax matrices in the most general case are 

2L = -\{U' l U) 2 + AU + BU 2 + k[U' 1 U, U] + h.c. (20) 

M = -\U~ l U + kU - h.c. (21) 

where h.c. denotes hermitian conjugate. Notice that L is hermitian while M is antihermi- 
tian. Relations (17) and (18) remain unchanged. 
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Alternatively, we could write the "right-Lax pair" matrices L, M, by substituting 
UU~ l for U~ l U in (20) and (21) and flipping the sign of the first term in M. The two 
choices are really equivalent, connected through a unitary transformation generated by U 
itself. 

The above are true without any assumptions for the constant commutator P. Assum- 
ing, now, that P is of the form (4) , we can express the conserved quantities I n in terms of 
particle coordinates x n and momenta p n = x n . To see this, perform a (time-dependent) 
unitary transformation on U which brings it to a diagonal form while it doesn't change I n . 
Then, using (4) , we see that the matrices of the system take the form 

Uij = <);,<'■'■* , Pij = a{8ij - 1) 



i - ia(l — 5jn) 

(g-'ioa = Hi/in + gg^n (22) 



The J n , therefore, can be expressed using (22) in terms of particle quantities and constitute 
N independent integrals of motion for the system of particles in the potential (8) interacting 
through inverse sine square two-body potentials of strength a 2 . I\ is the total energy of 
the system. I2 has the form 



(23) 



where we defined 



and used the identity 



ay = 2 sin - — 3 , B = \B\e l P (24) 



E — ; = ( 25 ) 

{i,j,k,l distinct} J 



We shall not give here the explicit form of the rest of I n in terms of x^ pi, which is quite 
complicated. Their independence is obvious from the fact that, just as in the hermitian 
case, I n is a polynomial in pi of order 2n with highest order term of the form J2iPf n ^ an d 
such terms cannot be obtained from lower order terms for n < N. 
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The above Lax pair matrices go over to the Lax matrices of the hermitian problem in 
the limit of infinite radius, under the simultaneous scaling (9) and (11) for the coordinates 
and the potential, and the rescaling 



In the case c = 0, though, there is no scaling limit, expressing the fact that the purely 
cubic potential does not admit a Lax pair formulation although it is still integrable. 

It should also be noted that, just as in the hermitian case, there are two distinct Lax 
pairs, corresponding to the two possible choices of sign for k in (18). Again, this has no 
impact on the conserved quantities I n since the two choices give essentially the same set 
of conserved quantities, modulo traces of P n , which are trivial nondynamical constants 
in the space of eigenvalues. Note also that, in the hermitian case, when the coefficient of 
the quartic term became negative the Lax matrix L became non-hermitian, indicating the 
runaway nature of the system. The integrals of motion, nevertheless, remained real. In 
the present case no such thing happens, since there can be no runaway behavior on the 
circle, and the Lax matrix (20) remains hermitian for all values of the coefficients of the 
potential. 

The special case B = is interesting: it corresponds to the quadratic hermitian model 
(i.e., the Calogero system) which is obtained as the large R scaling limit in this case. In 
particular, L = H just as in the hermitian case. Contrary to the Calogero system, however, 
where H is a constant matrix, in this case there is still some nontrivial unitary rotation 
with time, due to the nonzero curvature of the space. Therefore, the gauge matrix of the 
Lax pair M = —^U~^U, generating this rotation, does not vanish. At the scaling limit, 
of course, M goes to zero. 

It remains to show the fact that the above quantities are in involution, that is, that 
their Poisson brackets vanish. Again, as in the hermitian case, trying to use the explicit 
expressions of I n in terms of particle phase space variables and using 



L 




M 




(26) 



{ x ii Pj } — $ij 



(27) 



is completely hopeless. Instead, we shall work with the canonical structure of the original 
matrix problem and make use of the fact that the projection from the full matrix phase 
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space of U to the phase space of its eigenvalues is a hamiltonian reduction [6,7]. To see 
this, define from (1) a canonical momentum Pjj 

P u = ^ = -U^Ulf- 1 (28) 
U 8U 

Then the symplectic two-form u> is 

u = tr (dPu dU) (29) 

giving rise to canonical Poisson brackets for U and Pjj. For convenience, we shall work 
with the antihermitian matrix 

II = -P V U = U~ l U (30) 

and the symplectic one-form A, in terms of which we have 

A = -tr (n U~ l dU) , u = dA (31) 

Decompose, now, U and IT in terms of diagonal and angular degrees of freedom, namely 

U=V~ 1 AV, JI = V~ 1 (n + A~ 1 [A,A])V (32) 

with A and N diagonal unitary and antihermitian matrices, respectively: 

Ay = e Xi Sij , Nij = ip t 6ij (33) 

and A an off-diagonal antihermitian matrix which will be identified with the off-diagonal 
part of VV^ 1 by the equations of motion. Then, in this parametrization (31) becomes 

A = tr (-JVA _1 dA + (A^A^ 1 + A _1 AA - 2 A) dVV' 1 ^ 

= Y,Pi dx * + tr (iPV~ X dV^ ^ 

i 

So A decomposes into two non-mixing parts, the first one being the particle phase space 
one-form leading to (27) and the second part identifying P as the variable canonically 
conjugate to the angular degrees of freedom. Therefore, the Poisson brackets of quantities 
expressible in terms of and pi can equally well be evaluated using the full matrix phase 
space Poisson structure. 

To show the involution of I n it is much more convenient to use the geometric formu- 
lation of symplectic manifolds, rather than work explicitly with Poisson brackets. In this 
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approach, we map to each function / on the phase space a vector field vj on the phase 
space through the relation 

<Vf,u>=df (35) 

where < , > denotes internal product (contraction). The matrix differential operators Sji 
and 8jjU form a basis for the vector fields, satisfying 

< (S u )ij, dU kl >=< (SuU)ij, (U- 1 dU) kl >= 5 u S jk , else <.,.>= (36) 

In this basis, we can express the vector field Vf as 

v f = tr^^n + v 2 S u U) (37) 

Then the Poisson bracket of any two functions /, g can be calculated as 

U,9} = M9) = ^m + v2 W u ) (38) 

Therefore, defining v n = vj n , we have 

{I n ,Im} = Vn(Im) = mti (L m - l v n {L)) (39) 
On the other hand, v n satisfies 

< Un> oj >= dl n = t?(L n - l dL) (40) 
By expressing L in (20) in terms of IT and U and expanding dL as 

dL = ^2R l dHS l + s ^R l U- l dU S l (41) 

i i 
and using the expression for to derived from (31), that is 

uj = tr{-dW l dU + Il{U' l dU) 2 ^ (42) 

we get from (40) 

v n = -n^tr(^L n - 1 J R 2 % + [S i L n ~ 1 R^I^Sjj - S^R^ll) (43) 

i 

Finally, substituting (43) in (39) we obtain 

{I n ,I m } = -nm^tifl^^RiSjL^RjSi + L^R.SjL^RjUS^ - (m # n) (44) 

hi 
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The rest is a tedious algebraic exercise, calculating the explicit expressions for R^, Si, Rj 
and Sj and evaluating (44) . We shall omit all detail and simply give the result 

{J m , /„} = mn(n 2 - B) tr (l™- 1 !! L n ~ l \Sl, U]} + h.c. - (m ^ n) (45) 

Choosing then k as in (18) the Poisson brackets (45) vanish and the I n are in involution. 
We have therefore proved the full integrability of the system. 

The expression for L can be brought to a particularly suggestive form in a somewhat 
special case. Define the matrices 

A = U + Y(U) , A* = -U + Y(U) , with Y(U) = kU + k*U ^ + A , A = real (46) 

Then L = ^A^A is exactly the previous Lax matrix (20), with B and k properly connected. 
The coefficients of the potential, however, are constrained to obey the relation 

^2 = real ( 47 ) 

The situation is again analogous to the hermitian case, where such an expression for L 
also existed for a similarly constrained potential. Note that, in this case 

/„(«, A) = tr(l A^A) n = tr(l AA^) = /„(-«, -A) (48) 

which shows that the two possible choices of sign for k in (18) give exactly the same 
conserved quantities. The proof of involution of the I n is also simplified in this case. 

The above results can readily be extended to the case of positive definite matrices X. 
Such matrices can be expressed as 

X = e M , M hermitian. (49) 

They differ therefore from unitary matrices by a mere analytic continuation, and can 
be thought as unitary models with imaginary radius R. We can explicitly perform this 
continuation in the expressions for the Lax pair matrices, by putting R = i in (9), (10), (11) 
and (26), which amounts to the redefinitions 

U ^ X , -L 

a\ ± ia,2 — > -(oi ± CL2) , h ± ib^ — > -{b\ ± 62) (50) 
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and now k and k* are distinct, satisfying 

n 2 = -b\ - b 2 , k* 2 = ~h + b 2 



(51) 



The new potential of the system becomes 

V(x) = a\ cosh x + a 2 sinh x + b\ cosh 2x + 62 smri % x (52) 

Notice that, for (52) to be bounded from below, we must have 

bi ± b 2 > (53) 

and for such values k and k* both become imaginary. The Lax matrix becomes, explicitly 

L = \ll 2 + W(X) + ik 1 [U,X] + iK 2 [U,X- 1 } (54) 

where we put k = 2in\, k* = 2iK 2 - The momentum IT = X~^X now is not antihermitian 
any more, but rather satisfies 

nt = XliX 1 (55) 
Therefore, although L is not hermitian, it satisfies 

L f = XLX' 1 (56) 

which ensures that the conserved quantities I n are real. Alternatively, we can directly 
make the analytic continuation 

x % -> ixi , v% -»• m (57) 

into the expressions for I n and obtain the new expressions. These will be conserved quan- 
tities for particles on the line, in an external potential of the type (52) interacting through 
two-body potentials of the type 

Systems of particles with the above two-body potential and no external potentials were 
known to be integrable. As we proved, they remain integrable in the presence of potentials 
of the form (52). 
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In conclusion, we see that the above unitary system (and its positive definite ana- 
lytic continuation) is integrable and closely parallels the hermitian case. It is interesting 
that integrability again stops at potentials with only their two lowest nontrivial Fourier 
coefficients nonvanishing and does not extend to arbitrary periodic potentials. Obvious 
attempts to generalize the result to higher potentials, e.g., through a construction similar 
to (46) with arbitrary hermitian Y(U), fail. This may be just a shortcoming of the method 
used. If, on the other hand, it turns out that integrability indeed stops at this level, it 
would be interesting to understand the deeper reason for that and the special significance 
of these potentials. 

We should stress that, just as in the hermitian case, integrability does not necessarily 
imply complete solvability. The significance of the results is that the problem in the 
presence of the interparticle interaction remains as solvable as in the case of decoupled 
particles. The solution of these systems could conceivably be achieved in terms of the 
solutions of the one-body problem (which in general involves elliptic functions). The proof 
of integrability in the quantum domain is also of interest, especially in view of the conjecture 
put forth in [13] that the inverse square two-body potential on the line and the inverse sine 
square potential on the circle simply endow the particles with fractional statistics. These, 
as well as possible generalizations for other two-body potentials (e.g., of the Weierstrass 
form) remain interesting topics for further work. 
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